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Curves obtained with an analog computer are presented for the magnification factor 
versus frequency ratio of second-order systems with combined coulomb and viscous damping. 
The ranges of the parameters are as follows: Viscous damping ratio from 0.05 to 5.0, in 15 
steps; coulomb damping ratio from to 0.9, in 11 steps; and frequency ratio from to 2.0. 
Boundaries between regions with 0, 1, and 2 stops per half-cycle are also shown. 



1. Introduction 

The increasing importance of vibration measure- 
ments for both military and nonmilitary applica- 
tions has stimulated the production of a large variety 
of vibration -measuring instruments in recent years. 
This activity in design and manufacture has not 
been matched by progress in the analysis of the 
response of vibration instruments, or in the exten- 
sion or modernization of previous analyses of their 
response. It is still frequently assumed, in the 
application of these instruments, that the response 
of an actual vibration instrument is identical with 
the response of an ideal instrument in which the 
damping is entirely "viscous," i. e., proportional to 
the velocity of the mass or "seismic element." The 
actual response may differ significantly from the 
ideal response if the damping varies with other than 
the first power of the velocity, or is dependent on 
displacement, or if any appreciable amount of 
coulomb damping is present. This last type of 
damping, sometimes referred to as dry or sliding 
friction, exerts on the moving element a force that 
is constant in magnitude and always acts in the 
direction opposite to the velocity of the moving 
element. 

Coulomb damping is usually present in a mechani- 
cal second-order system such as is illustrated in 
figure 1. Here, in addition to the coulomb damping, 
a spring force and a viscous damping force are shown 
as acting upon the seismic element. The motion of 
such a system has been analyzed by Den Hartog 3 
in a 1931 paper in which results are presented 
graphically in a series of figures. For application to 
vib ration i nstruments in present use, Den Hartog's 

1 This work was conducted under a program of basic instrumentation research 
and development sponsored by the National Bureau of Standards, the Office of 
Naval Research, the Air Research and Development Command and the Atomic 
Energy Commission. 

2 Now at the Missile Systems Division .Lockheed Aircraft Corp., Van Nuys, 
Calif. 

•J. P. Den Hartog, Forced vibrations with combined coulomb and viscou 
friction, Trans. Am. Soc. Mech. Eng. 53, APM 53-9, 107-115 (1931). 



results require extension to larger values of viscous 
friction and a larger range of frequencies. 

The present paper extends these results through 
the ranges applicable to such inertial instruments as 
accelerometers and jerkmeters. These inertial in- 
struments are usually designed for viscous damping 
between 0.6 and 0.7 of critical damping, and are 
usually employed between zero frequency and a 
maximum frequency less than the resonance fre- 
quency (of the mass and spring) . Particular empha- 
sis is therefore placed in this paper upon results for 
viscous damping between 0.6 and 0.7 of critical 
damping, and for frequencies between zero and the 
resonance frequency. 

Results are presented graphically in figures 2 to 
17. These curves extend from zero to twice resonance 
frequency, and are presented so as to be readily 
applicable to vibration instruments used below their 
resonance frequency. They were determined without 
particular difficulty through the use of an analog 
computer. The original program called for the 
computation of an additional set of curves for 
frequencies between twice the resonance frequency 
and 50 times the resonance frequency. These 
curves would have been applicable to devices such 
as seismometers and velocity meters, which are 
ordinarily employed at frequencies above their 
resonance frequency. Unfortunately, the analog 
computer accuracy proved poor at these high fre- 
quencies. The computation of displacement at high 
frequencies consisted essentially of the double inte- 
gration of the difference of a sinusoidal forcing func- 
tion and the coulomb friction. When the analog 
computer attempted this double integration, its results 
contained appreciable errors caused by drifting and 
hunting. The size of these errors caused a termina- 
tion of the high-frequency computations. Hence, 
only curves applicable to vibration instruments 
usually employed below their resonance frequency 
are presented in this paper. 
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2. Differential Equation for Vibration 
Instruments 

figure 1 shows a typical vibration instrument 
with its frame rigidly attached to a sinusoidally 
vibrating structure having a motion X cos ut. The 
displacement, x, of the seismic element relative to the 
frame of the instrument satisfies the differential 
equation : 



Mx+Cx+Kx±F=MXa> 2 cos <rt, 



(1) 



where (see fig. 1) 



M=mass of seismic element. 
C= viscous damping force per unit velocity. 
K= restoring force exerted by the spring per 

unit displacement. 
F= coulomb damping force, assumed constant 
in magnitude and changing sign so as 
always to oppose the motion. 
X= amplitude of forcing motion. 
co=27r/=angular frequency of forcing motion. 

A number of equations in Den Hartog's paper give 
a complete literal solution of this equation with the 
right-hand side of the equation replaced by a term 
P cos (co£+<p), where P is the amplitude of a "periodic 
disturbing force. " For practical application Den 
Hartog expressed his results by a set of computed 
curves of magnification factor versus frequency ratio 
for seven values of viscous damping ratio between 
and 0.5 of critical damping, and for a number of 
values between and 0.9 of the ratio of coulomb 
damping force to peak disturbing force. Most of 
Den Hartog's curves extend over a frequency range 
of 0.5 to 2.0; the one for zero viscous damping ex- 
tends over a frequency range of to 2.0. Curves for 
zero coulomb damping (i. e., viscous damping only) 
are also available elsewhere. 4 In the present paper, 
all the curves extend over the frequency ratio range 
between and 2.0. They are given for 15 values of 
viscous damping ratio between 0.05 and 5.0, and for 
1 1 values of coulomb damping ratio between and 
0.9. 

3. Summary of Results 

The curves shown in figures 2 to 17 are drawn 
through points plotted by an analog computer. Each 
figure presents a family of curves for a constant value 
of C/C c , the ratio of vis cous damping C, to critical 
viscous damping, C C =2^KM. Each curve of the 
family is a plot of magnification factor or acceleration 
response x ool/X versus the frequency ratio, a>/co w , 
where x is t he maximum value of the displacement x, 
co n ='\lK/M is the undamped resonance frequency of 



^See, for example, C. S. Draper, W. McKay, and S. Lees, Instrument engineer- 
ing. II. Mathematics (McGraw-Hill Book Co., Inc., New York, N. Y., 1953). 



the mass-spring system, and X=o) 2 x is the peak 
acceleration of the forcing motion. For inertial 
instruments, the term "acceleration response" seems 

appropriate for x o) n 2 JX because this quantity may be 
interpreted as the product of two important instru- 
ment characteristics: Xq/X, the response (displace- 
ment) of the seismic element per unit acceleration 
applied to the instrument, and c4, the resonance 
frequency of the instrument, to which instrument 
sensitivity is inversely proportional. 

Each curve of the family is drawn for a particular 

value of the parameter a F /X, where a F is the mini- 
mum acceleration that must be applied to the instru- 
ment to overcome coulomb friction and produce 

motion. This parameter a F /X is equivalent to the 
parameter F/P used by Den Hartog. Either may be 
obtained from the other by employing the relations 

F=Ma F and P=MX. 

The dashed line plotted in each figure is a boundary 
line between steady-state motions with no stop and 
one stop per half-cycle. In the region above this 
boundary, the steady-state motion is continuous 
without stops. In the region immediately below this 
boundary, the motion stops once each half-cycle. 

Only these two regions, nonstop and one-stop 
motion, were computed for C/C c > 0.2 and u/w n >0.1. 
For C/C c < 0.1, four regions characterized by no stop, 
one stop, two stops, or more than two stops each 
half-cycle were computed. When more than two 
regions occur, the magnification factor versus fre- 
quency ratio curves become very irregular for co/co w 
less than about 0.4. For w/w w <0.1, larger and larger 
amounts of computing time would be required for 
accurately determining the magnification factor as 
O/Ce becomes smaller and smaller. Hence for 
C/C c =0.2, 0.1, and 0.05, the curves are terminated 
at a?/co n =0.1. 

4. Accuracy of Results 

The curves of figures 2 to 17 are plotted to permit 
reading the magnification factor to within 1 percent 
of the full scale of each figure. These curves were 
drawn through points plotted by an analog computer. 
Over-all accuracy of the points plotted by the analog 
computer, which consisted of elements with an ac- 
curacy of ±0.1 percent of rated output, is estimated 
at ±0.3 percent of the full-scale magnification factor. 
A number of check points computed manually and 
by the computer agreed with each other to within 
±1.0 percent. Total computer errors and drafting 
errors are believed to introduce, in any figure, inac- 
curacies of no more than 1 percent of the full-scale 
magnification factor of each figure. 

It is believed that greater accuracy in computa- 
tion would have resulted from the use of a digital 
computer that employed Den Hartog's literal equa- 
tions. However, an estimate of total costs showed 
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analog computation to be more economical. One 
reason was the availability in the analog computer's 
auxiliary equipment of a plotting board, which could 
partially process the results of the computer study by 
plotting points of the magnification-factor curves. 
The other important important advantage of the 
analog computer was its use of eq (1) rather than the 
involved literal equations of Den Hartog's paper, 
which would have been employed by the digital com- 
puter. Use of the analog computer avoided the ex- 
penditures necessary for a careful scrutiny of these 
equations for possible misprints. 

5. General Effects of Coulomb Friction Upon 
Instrument Response 

The presence of coulomb friction makes instru- 
ment response dependent on the amplitude of the 
sinusoidal excitation. This amplitude dependence 
is shown graphically in figure 18, which is a cross-plot 
of the data for C/C c = 0.65, as given in figure 10. In 
figure 18, magnification factor is plotted versus 

Xja F for various constant values of co/o>„. 

The asymptotes shown as horizontal lines in the 
right-hand portion of figure 18 represent the magni- 
fication factors at vibration amplitudes sufficiently 
large to make negligible the effect of coulomb friction. 
An inspection of the figure shows that, for C/C C =0.Q5, 

X\a F must be >25 if the actual response is to be 
within 5 percent of the asymptote for all values of 

0)<OJ n . 

Rules similar to this can be stated for other values 
of C/C c . For CIC C >0.3 and for <)<«/«„< 1.0, a 
general rule of thumb for estimating the error intro- 
duced by coulomb friction is to assume that this 
error is equal to the error at co/co n =0. According to 

this general rule, a given percentage value of a F /X 
results in that same percentage reduction of magnifi- 
cation factor from its ideal (a F /X=0) value. This 
rule is a good approximation for the reduction in re- 
sponse for C/O c in the range 0.6 to 0.7. As the value 
of C/C c moves further and further outside this range, 
the approximation becomes poorer and poorer. 
However, the damping in most inertial instruments 
is between 0.6 and 0.7 of critical damping; hence, 
this rule is useful for practical inertial instruments. 

A physical example of the effect of coulomb frac- 
tion upon instrument response is shown in figure 19. 
This figure reproduces an oscillogram from two 
different velocity meters mounted back to back and 
subjected to the same sinusoidal forcing function. 
The upper trace shows that one instrument is stop- 
ping once each half-cycle as the result of coulomb 
friction. The lower trace shows the other to be 
exhibiting essentially ideal behavior. 

Figure 20 shows the time variation of displacement 
and velocity for various values of a F /X, for a low 



value of viscous friction (C/C c =0.1), and a low value 
of frequency ratio (co/co n =0.1). These records were 
made during the computer study. Each set of two 

traces is for a different value of a F /X. At a F /X=0.3, 
one-stop motion occurs, as evidenced by the long dwell 
of the velocity trace on the (central) zero line. For 

a F /X equal to 0.5, the second trough of the damped 
sinusoidal oscillation in the velocity trace approaches 

the zero line. At a F /X= 0.525, this trough reaches 
the zero, and the motion stops twice during each 

half-cycle. Further increase of a F /X to 0.55 fur- 
nishes an additional restraint, which prevents 
another start after the motion stops at nearly the 
same point of the velocity trough. Thus we have 
again one stop per half-cycle. The motion remains 
one-stop until displacement and velocity decrease 
sufficiently for a stop to occur at the first trough of 
the velocity, giving two-stop motion once more. 
This occurs at approximately a F /X=0.8. However, 

an increase of a F /X to 0.9 furnishes the additional 
restraint necessary to prevent another start. Thus 
motion again becomes one-stop. 

AsC/Cc decreases from 0.1, the damped sinusoid of 
velocity becomes less damped and it is conceivable 
that for co/w n <0.1, motion may stop not only at the 
third, second, or first trough of the velocity, but also 
at a fourth, fifth, or higher-order trough of the 
velocity. In such cases, the response of the instru- 
ment becomes very irregular. 



The analog computer study was performed by W. 
McCool and B. Zimmerman of the Naval Research 
Laboratory. The authors thank Messrs. McCool 
and Zimmerman for their suggestions and excellent 
cooperation in this study. M. Abramowitz and R. 
Dressier supplied information on the use of the digital 
computer. W. A. Wildhack continually furnished 
stimulating ideas during the investigation. A. Bucek 
prepared the curves and illustrations shown in this 
report. 
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Figure 1. Typical mechanical vibration instrument. 
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Figure 2. Instrument response for C/C e =0.05. 
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Figure 3. Instrument response for C/C c =0.05. 
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Figure 4. Instrument response for C/C c =9.10. 
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Figure 5. Instrument response for C/C c —0.gO. 



il 




0.6 0.8 1.0 1.2 1.4 

Frequency Ratio ^A^r\ 



2.0 



Figuke 6. Instrument response for C/C c = 0.30. 
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Figure 7. Instrument response for C/C c = 0.40. 
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Figure 8. Instrument response for C/C c =0.50. 
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Figure 9. Instrument response for C/C c =0.60. 
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Figure 10. Instrument response for C/C c =0.6b~. 
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Figure 11. Instrument response for C/C c =0.70. 
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Figure 12. Instrument response for C/C C =0.7S. 
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Figure 13. Instrument response for C/C c =0.80. 
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Figure 14. Instrument response for C/C c =0.90. 
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Figure 15. Instrument response for C/C c =1.0. 
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Figure 16. Instrument response for C/C c =2.0. 
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Figure 17. Instrument res-ponse for C/C c =5.0. 
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Figure 18. Effect of amplitude upon instrument response. 
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Figure 19. Experimental response of two different 
vibration instruments. 
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Figure 20. Computed response of vibration instrument. 

Washington, April 3, 1956. 
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